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2 ■ Lyman-alpha (Lya) is one of the dominant tools used to probe the star-forming 

' galaxy population at high-redshift (z). However, astrophysical interpretations of data 

^ I drawn from Lya alone hinge on the Lya escape fraction which, due to the complex 

radiative transport, may vary greatly. Here we map the Lya emission from the local 
luminous blue compact galaxy Haro 11, a known emitter of Lya and the only known 
candidate for low-z Lyman continuum emission (LyC). To aid in the interpretation 
' we perform a detailed UV and optical multi-wavelength analysis and model the stellar 

I population, dust distribution, ionising photon budget, and star-cluster population. We 

, use archival X-ray observations to further constrain properties of the starburst and 

^SJ ' estimate the neutral hydrogen column density. 

The Lya morphology is found to be largely symmetric around a single young star 
forming knot and is strongly decoupled from other wavelengths. From general surface 
' photometry, only very slight correlation is found between Lya and Ha, E{B — V), 

, and the age of the stellar population. Only around the central Lya-bright cluster do 

we find the Lya/Ha ratio at values predicted by recombination theory. The total Lya 
escape fraction is found to be just 3%. We compute that ^ 90% of the Lya photons 
that escape do so after undergoing multiple resonance scattering events, masking their 
' point of origin. This leads to a largely symmetric distribution and, by increasing 

I the distance that photons must travel to escape, decreases the escape probability 

significantly. While dust must ultimately be responsible for the destruction of Lya, it 
plays little role in governing the observed morphology, which is regulated more by ISM 
kinematics and geometry. We find tentative evidence for local Lya equivalent width 
in the immediate vicinity of star-clusters being a function of cluster age, consistent 
with hydrodynamic studies. We estimate the intrinsic production of ionising photons 
and put further constraints of ~ 9% on the escaping fraction of photons at 900A. 

Key words: galaxies: evolution - galaxies: star clusters - ultraviolet: galaxies - 
galaxies: individual: Haro 11 
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1 INTRODUCTION 

Four decades ago iPartridge fc PeeblesI (|l967l ) discussed the 
prospects of identifying 'primeval' galaxies (i.e. galaxies 
forming their first generation of stars), using both the Ly- 
man decrement and Lyman-alpha (Lya) emission line as 
observational tracers. Both methods rely upon the galaxies 
hosting numerous young, massive stars producing a strong 
radiation field in the far ultraviolet (UV). The absorption 
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of this radiation field bluewards of tlie Lyman absorption 
edge results in the Lyman-break phenomenon, while the re- 
processing of the absorbed photons in astrophysical nebulae 
results in the superimposition of strong hydrogen recom- 
bination lines on the galaxy spectra. While it's likely that 
both features are present in the spectra of high-redshift (z) 
starbursts, from a survey perspective they compete. Lyman- 
break candidates are identified in multi-broadband imaging 
surveys where large ranges in redshift (and therefore cosmic 
volume) can be probed, but are biased towards the detec- 
tion of galaxies with FUV continua towards the brighter 
end of the luminosity function (LF) - less numerous in the 
hierarchical galaxy formation scenario. On the other hand, 
because emission lines concentrate a large amount of energy 
in a very small spectral region, sources with much fainter 
continuum can be uncovered. Unfortunately, isolation of a 
line requires either a spectroscopic surveys which typically 
probe narrow fields of view because of the slit dramatically 
reducing the size in one dimension, or narrowband imaging 
which probes only a small range in redshift. The advent of 
large diameter reflectors, efficient optics, and sensitive pho- 
tometric devices has resulted in both techniques enjoying 
much success in recent years. 

Young galaxies with low dust contents can be ex- 
pected to produce Lya emission with high equivalent widths 
{Wi^ya), ~ 240A for solar-like meta llicities (Z) and stan - 
dard initial mass functions (IMF) jCharlot fc Falll Il993h . 
This value is predicted to increase significa ntly as metallic - 
ities approach the population ill domain (ISchaereij l2003h . 
Since the Lya rest wavelength lies in the FUV, the line 
is still observable from the ground in the optical domain 
at z ~ 6 and beyond making it, in principle, the ideal 
tool by which to identify star-for ming galaxies in the early 
universe. Despite a rocky start (|Pritchetlll993l ) Lya emit- 
ters (LAEs) are showing up en masse in high-z surveys 
and perhaps their obs erved and predicted nu mber-densities 
are now converging (|Le Delliou et al ] 120061) . Lya IS now 
being used to p ut constraints on the final stages of cos- 
mic r eionisation (iMalhotra fc Rhoadsl |2003 : iDiikstra et al.l 
I2OO6I ) and explor e high-redshift large scale structure and 
galaxy clustering llHamana et al ]|2004l : iMalhotra et~al]|2005l : 



iMuravama et ah l2007h . In addition, numerous sources have 
been detected t hrough narrowband imaging techniques with 
very high lyr.v^fMalho tra &: Rhoadsll2002l : IShimasaku et al.l 
I2OO6I ). perhaps indicative of top heavy IMFs or extreme 
stars, making such targets ideal to search for signatures 
of the so-far illusive population ill stars. However, so far 
spectroscopic observations have failed to detect the strong 
HeiiA1640A feature expected from pop iii objects (eg. 
iDawson et~aLll2004 [Nagao et al.ll2007l '). Finally, Lya is also 
being used to estimate the cosmic star format i on rate density 
at th e highest redshifts l|Fuiita et al.l l2003l : lYamada et al.l 
I2OO5I ). Such star formation rates (SFR ) are typically esti - 
mated assuming case B recomb ination (|Brocklehurstlll97"ll ) 
and the Ha SFR calibration of iKennicuttI 1 19981 ). although 
the large spread in SFR(Lya) vs. SF R(FUV) and consistent 
underestimates from SFR(Lya) (eg. iMuravama et al1l2007l ) 
suggest that such a calibration should be used with caution. 
This caution must be extended to all cosmological studies in 
light of the fact that only a fraction of UV-s elected high- 2: 
targe ts show a Lya feature in emission (e.g. IShaplev et al.l 
l2003h . 



The complexities of using Lya as a cosmological tool 
result from the fact that it is a resonant line, and its po- 
tential cosmological importance has motivated a number 
of studies of star-forming galaxies at low-z and theoreti- 
cal models of resonant line radiative transfer. Early stud- 
ies with the International Ultraviolet Explorer (lUE) ini- 
tially mirrored the first results at high-z: Lya was typically 
weak or absent in local starbursts. This systematic weak- 
ening of Lya was fi rst attributed to dust absorption (eg. 
ICharlot fc Faiilll99ll ) and an anticorrelation between VFlvq 
and metallicity was found in the /[/i? sample (I Chariot fc Falll 
|T9_93;). However, the damped Lya absorption seen in some 
local starbursts, and the failure of dust corrections to recon- 
cile Lyg with the fiuxes p redicted by recombination the- 
ory l|Giavalisco et al.1 Il996l ). is indicative of selective at- 
tenuation of Lya. This is clearly exemplified by the fact 
that the most metal-poor galaxies known at low-z (iZw 18 
and SBS 03350-52) are shown to be damp ed Lya absorbers 
(|Kunth et al.lll994l : iThuan fc Izotovlll997l ) in spectroscopic 
observations from the Goddard High Resolution Spectrograph 
(GHRS) and Space Telescope Imaging Spectrograph (STIS). 
This is to be expected if the starburst is enshrouded by a 
static layer of neutral hydrogen that is able to resonantly 
trap Lya, thereby greatly increasing the path-lengths of 
the photons in order to escape the host, and exponentially 
increasing the chance of their destruction by dust grains 
(Neufcld . 19901 Further G/fTiS observations (iKunth et al.1 
Il998h showed that when Lya is seen in emission, it almost 
systematically shows a P Cygni profile and systematic veloc- 
ity offset from low ionisation state metal absorption features 
in the neutral ISM, suggesting that an outfiowing medium 
is an essential ingredient in the formation of the line profile 
and Lya e scape physics. Similar results have also be en found 
at high-z l|Shaplev et al.ll2003l : pTapken et al.ll200'it). Hydro- 
dyna mic models of expanding bubbles ( Tenorio-Tagle et al.l 
Il999h predicted Lya line profiles to follow an evolutionary 
sequence starting with pure absorption at the earliest times, 
developing through a pure emission phase into P Cygni pro- 
files as the ISM is driven out by mechanical feedback from 
the starburst, and fading back i nto ab sorption at late times. 
This allowed iMas-Hesse et al.1 l|2003l ) to reconcile the vari- 
ety of Lya profiles observed at low-z in suc h an evolutionary 
sequence. Further theo retical studies (e.g. lAhn et al.|[2003l : 
IVerhamme et all I2OO6I ) have shown how a wide variety of 
P Cygni-like and asymmetric profiles can develop, depending 
on the ISM properties and geometry. The line form ation be- 
come s more complex still if the ISM is multiphase (|Neufeldl 
I1991I : n ansen fc OhI I2OO6I ) and certain physical and kine- 
matic configurations may lead to increases in the fraction 
of escaping Lya photons compared to non-resonant radi- 
ation, thereby effectively boosting Whya- Recent advances 
have been made in the understanding of observed line pro- 
files with the implementation of full 3-D codes with arbitrary 
distribu tions of gas, ionisation, temperature, dust, and kine- 
matics (jVerhamme et aLll2006l ). 

Since nebular ionisation does not occur in situ with ion- 
ising sources, recombination line imaging (e.g. Ha) may re- 
veal a morphology different from that found by imaging of 
the stellar continuum. This phenomenon may be much more 
significant for Lya photons which resonantly scatter. While 
Lya and Ha have the same points of origin, it is likely that 
internal Lya scattering events may cause Lya photons to 
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be emitted far from their production sites, and not be spa- 
tially correlated with Ha or other non-resonant recombina- 
tion lines. Targeted spectroscopic observations are therefore 
liable to miss a fraction of the emission and, while contain- 
ing neither frequency nor kinematic information, Lya imag- 
ing becomes an invaluable complement to the spectroscopic 
studies. This was the motivation for our imaging survey 
of local starbursts us ing th e Advanced Camer a for Surveys 
(ACS ) l|Kunthet al.ll2003l : iHaves et al.1l2005l : lOsthn et al.1 
I2OO7I in preparation). This is a truly unique dataset for a 
number of reasons. Firstly the angular resolution of the ACS 
allows us to map the Lya morphology on scales of 5-15 pc; 
2-3 orders of magnitude better than typical studies at high- 
z. The addition of Ha (absent in almost all high-2 studies 
due to an inconvenient rest wavelength) allows us to quan- 
titatively examine the decoupling of Lya from non-resonant 
lines and estimate the global escape fractions. Multiband 
UV and optical data allow us to map dust reddening, stel- 
lar ages and masses, ionising photon production, and other 
properties of the host, all at the same resolution as Lya. 
The first detail ed Lya imaging st udy of a local starburst, 
ESO338-IG04 l|Haves et al. 1120051 ') found emission and ab- 
sorption varying on very small scales in the central starburst 
regions, and little or no correlation with the FUV morphol- 
ogy. The starburst is surrounded with a large, diffuse, low 
surface brightness Lya halo that contributes ~ 70% to the 
global Lya luminosity, resulting from the resonant decou- 
phng and diffusion of Lya. The total escape fraction was 
found to be just ~ 5%, implying any global values (eg. SFR) 
that would be estimated from Lya alone would be seriously 
at fault. 

Feedback processes from star-formation are capable of 
driving galaxy-scale 'superwinds' that shock heat and accel- 
erate the ambient medium and circumnuclear gas, resulting 
in large-scale, diffuse X-ray nebulae (Hcckman ct al. 19901). 
Typically, hot X-ray emitting regions exhibit a tight mor- 
phological correlati on with with warm gas as probed by op- 
tical emission lines l|Grimes et ahlfiooH ). In deed, such X-ray 
nebulae are near-ubiquitous in starbursts ([Strickland et al] 
I2OO4I ') and have been calibrated as tracers of s tar-formation 
rates llRanalh et al.ll200l : IColbert et al]|2004l ) . iGrimes et al] 
(|2005l ) have also shown that observed X-ray spectra are well 
reproduced by a single thermal plasma over a range of galaxy 
classifications spanning dwarfs, discs, and ultra-luminous in- 
frared galaxies (ULIRG), although the central regions of 
some ULIRGs require an additional power-law component. 
Since the diffuse X-ray component can be so valuable for an 
understanding of the wind properties, a nd feedback is an es - 
sential ingredient in the escape of Lya (|Kunth et alj[l99^ . 
X-ray information provides a valuable supporting dataset 
for a detailed investigation of Lya. 

In this article we turn our attention to another tar- 
get in our sample, the well-known, lumin ous (Mb = -20.5), 
low m etallicity {\og{0/H) + 12 = 7.9; iBergvaU fc OsthnI 
|2002| ) blue compact galaxy (BCG). It exhibits a complex 
morphology consisting of three main star-for ming conden- 
sations and an unrelaxed kinematic structure (|6stlin et al.l 
[l999. . .2001 ). suggestive of a dwarf galaxy merger. It is ac- 
tively star-for ming, exhibits a large number of bright young 
star clusters (lOsthnI |2000D . is a known emitter of Lya 
jKunth et al.l [l99^ ). and the only k nown local candidate 
emitter of Lyman continuum (LyC) l|Bergvall et al.l I2OO6I : 
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iGrimes et ai]|2007l ). The FUV continuum luminosity puts 
Haro 11 at the brighter end of the di stribution of Lya emit- 
ters at z ~ 3.1 (|Gronwall et al.ll2007l ). We utilise HST/ACS 
images in the FUV {Solar Blind Channel; SBC), to examine 
Lya and nearby continuum, and broadband images in the 
UV and optical, and narrowband Ha {High Resolution Cam- 
era; HRC, and Wide Field Camera; WFC) to examine dust, 
ages in the stellar population, the star cluster population 
as a whole, and estimate the LyC production. We use deep 
ground-based narrowband images in Ha and H/3 in order 
to estimate extinction in the gas phase. X-ray observations 
of Haro 11 have be en obtained using the Chandra satellite 
IGrimes et alJlioOTl ) but currently their status is still propri- 
etary. We hence also exploit serendipitous off-axis observa- 
tions from the Chandra and XMM-Newton X-ray observato- 
ries to study the wind properties and internal photoelectric 
absorption. 

The article is arranged in the following manner: in Sec- 
tion [2] we describe the observations and data reductions; in 
Section[3]we present the results; in Section|3]we analyse dis- 
cuss the results; and in Section[S]we present our concluding 
remarks. We assume a cosmology of Ho = 72 km s~^ Mpc^^, 
SIm = 0.3 and Qa = 0.7 throughout. The redshift of Haro 11 
is taken to be 0.020598 from NED, corresponding to a lumi- 
nosity distance of 87.1 Mpc. 

2 OBSERVATIONS, REDUCTIONS, AND 
DATA PROCESSING 

2.1 HST Ultraviolet and optical 

2.1.1 Observations 

This study makes use of images from all three channels of 
the ACS onboard HST: the SBC for the FUV; the HRC for 
the 22OOA and (7— band; and the more sensitive WFC for im- 
ages in the optical domain. The details of the observations 
and post-reduction processing of the HST UV and optical 
dataset are described in a companion paper, lOstlin et al-l 
(2007). Observations were performed in the bandpasses 
listed in Table [T] Briefiy, F122M and FI4OLP correspond 
to Lya on-line and continuum, respectively, while FR656N 
covers Ha on-line and F550M measures line-free continuum 
bluewards of Ha. The remainder are broadband continuum 
observations, selected in order to avoid the strongest emis- 
sion lines. 

All images are 'drizzled' using the MULTIDRIZZLE 
task in IRAF/STSDAS onto the same pixel sampling scale 
(0.025"/pixel) and position angle. The inverse variance 
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weight maps were saved from the drizzle process since they 
provide an estimate of the error in each pixeL Remaining 
cosmic rays, charge transfer tracks, and blemishes are re- 
moved from the CCD observations using the CREDIT task. 
Remaining band-to-band discrepancies in the astrometric 
ahgnment were rectified using the GEOMAP and GEOTRAN 
tasks. Since our dataset comprises observations between 
1200 and 9000A and utilises three different ^4 (75 channels, we 
address the issue of variations in the point spread function 
(PSF) of the images. PSF models of each band were gen- 
erated using the PSF task in DIGIPHOT/DAOPHOT, with the 
resulting models being used to convolve all images to the 
PSF of F550M (the broadest emission-line free PSF in our 
dataset). PSF models were re-computed for the convolved 
images and compared to that of F550M; all frames showed 
PSF full width at half maximum consistent with F550M at 
below the 5% level. 



2.1.2 Lya continuum subtraction and SED modeling 

As described in iHaves et al.1 (|2005l ). and in more detail in 
Hayes et al. (2007, in preparation), continuum subtraction 
of Lya using these filters requires more sophisticated tech- 
niques than most emission lines. Because the offline filter 
is removed from the online filter by AA/A = 0.22 and the 
FUV continuum evolves rapidly as a function of A, age, and 
E{B — V), it is imperative to understand the behaviour of 
the continuum between FI4OLP and F122M. Special care is 
also required because of the broad nature of the F122M fil- 
ter which has a rectangular width around 10% of the central 
wavelength, and transmits both the geocoronal Lya line and 
Milky Way Lya absorption profile. In addition, P Cygni pro- 
files may result in a reduction of the net emission due to the 
blue-side absorption cancelling some or a ll of the emission. 

We defined in iHaves et al. I (|2005l ) the Continuum 
Throughput Normalisation (CTN) factor as the factor that, 
for a given spectrum, scales the continuum flux sampled 
in FI4OLP to the flux that would be expected from con- 
tinuum processes alone in F122M. The procedure of es- 
timating CTN in each pixel utilises each of the images 
between F14-0LP and F8I4W, the filter throughput pro- 
files, and_S\AAng_Sta^ 

els llLeitherer et al.1 1 19991 : IVazguez fc Leithereij 120051 ). In 
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(|2005l ) we demonstrated that we could find 
non-degenerate best-fitting spectra if SED datapoints sam- 
ple the UV continuum slope and Balmer/4000A break. That 
way, for each pixel we fit burst age and E{B — V) using stan- 
dard minimisation. The method has been substantially 
developed and is described in detail in Hayes et al. (2007, in 
preparation) . We now use continuum-subtracted Ha to first 
map the contribution to the overall SED from nebular gas 
alone. Thanks to observations at V and I we are also able to 
constrain the contribution from any stellar population that 
underlies the current starburst. Essentially we measure and 
subtract the gas spectrum and fit age and mass in two stel- 
lar components, applying the same reddening for all three 
SEDs. For each pixel, we are able to find the age of the stel- 
lar populations and E{B — V), treating the nebular gas and 
two stellar components independently. With the best-fitting 
spectral reconstruction (composite Starburst99 spectra) at 
each HST pixel we then map the CTN factor and are able 
to reliably continuum subtract Lya. 



In Hayes et al. (2007, in preparation) we present ex- 
tensive simulations, designed to test the reliability of this 
methodology for an array of input spectra. We determine 
that the method must account for the presence of an under- 
lying stellar population and nebular continuum emission. 
Both these components may affect the U — B colour and 
result in bad fits and poor recovery of CTN in cases where 
they dominate the optical luminosity. Our method of recon- 
structing the nebular continuum and fitting multiple stellar 
components allows robust recovery of the FUV continuum 
features and CTN. An independent measure of the metallic- 
ity is also found to be a requirement although the metallicity 
of Haro 11 is well known. Details of the IMF are found not to 
have a detrimental impact when multiple stellar component 
fitting (i.e. a more complex star-formation history) is used. 
Overall we find that the software employed here is able to 
always recover input Lya equivalent widths to within 30% 
for 'weak' Lya emission (W^Lya = 10A) and to within 10% 
when the Lya line is stronger (WLya = 100A). With regard 
to the integrated fiuxe s and visual morphologies, we found 
in iHaves et al. I l|2005h that varying the parameters of the 
model spectra had a minimal effect on our results. Morpholo- 
gies were always indistinguishable by eye and integrated Lya 
fluxes self-consistent to within ~ 25%, even when pushing 
the parameters outside the regimes deemed reasonable for 
the galaxy in question. 

This procedure also provides numerous other outputs 
against which we can compare Lya. The full list of output 
maps is: continuum subtracted Lya and Ha, CT A'^-factor 
at Lya {F122M/F140LP), continuum flux densities at 900A 
(/900), I5OOA (/i50o), and 22OOA (/2200), E{B-V), the age 
of the two stellar components, the mass of the two stellar 
components, the stellar equivalent width in absorption of Ha 
and H/3, and ■ These allow us to compare Lya fluxes and 
equivalent widths with various local stellar ages, E{B — V), 
local star-formation rates, and estimate the ionising photon 
production. 

Haro 11 metallicity is known to be around 20% solar 
l|Bergvall fc 6stlinll2002l ) so for this study we adopt the evo- 
lutionary tracks generated with metallicity Z — 0.004. Since 
we are concerned with OB-dominated, young starbursts we 
use the tracks of the Geneva group. In the absence of a 
quantitative estimate, the IMF was assumed to follow that 
of Salpeter {a = -2.35) in the range 0.1 Mq to 120 Mq. 
Star-formation history was taken to be that of an instan- 
taneous starburst (single stellar population) with a mass 



2.1.3 Binning 

As described in the introduction, we can expect to see 
Lya features (either emission or absorption) around clus- 
ters where variations can be expected on small scales, and/or 
large-scale diffuse emission at low surface brightness. Typ- 
ically S/N per resolution element is high in the vicinity of 
strong continuum sources but signiflcantly lower (< 1) in 
the diffuse regions. Ordinarily, such problems are overcome 
by smoothing. While flxed kernel smoothing may improve 
S/N in diffuse regions it smears out spatial details in regions 
where S/N is good and smoothing would not be necessary. 
Adaptive smoothing may maintain some spatial resolution 
but does not necessarily conserve surface brightness. 
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Adaptive binning provides an alternative to the smooth- 
ing approach and f or this study we make u se of the Voronoi 
tessellation code of IPiehl fc Statleil ll20p6[). a ge neralisa tion 
of the Voronoi binning algorithm of Cappellari fc CopinI 
Voronoi tessellation overcomes the drawbacks of 
smoothing by binning together groups of pixels to conglom- 
erate resolution elements (frequently referred to as "spax- 
els"), recomputing the signal and noise in each new bin. Pix- 
els are continuously accreted until a threshold S/N has been 
met in each spaxel. The advantage of this adaptive binning 
method is that in high-S/A'^ regions, the spatial sampling 
remains high because spaxels are typically small, whereas 
in diffuse regions S/N is improved greatly. Diffuse emission 
regions, by definition, show variations over much larger spa- 
tial scales. Since each individual pixel is used exactly once, 
surface brightness is always conserved. 

The FI4OLP FUV science observation is assigned as 
the "reference" image, and used to generate the binning 
pattern, using the inverse variance map output by the 
STSDAS/MULTIDRIZZLE task to compute S/N. Pixels are ac- 
creted into spaxels until S/N = 5 has been met, with a 
maximum size of 40^pixels = ID" in order to speed up the 
process and prevent spaxels from growing arbitrarily large. 



2.1.4 Super star clusters 

Point-like sources are identified in the deepest of the opti- 
cal bandpasses, F435W and F550M using the DAOFIND task 
in IRAF. To be considered a detection, a cluster must be 
present in both of these bands. Each object was then in- 
spected by eye and any detections that were clearly spuri- 
ous were removed from the catalogue. For crowded fields 
PSF-photometry is preferred to standard aperture meth- 
ods in order to eliminate cross-contamination from neigh- 
bouring clusters. However, some clusters may be extended 
enough to be resolved by our observations, thus precluding 
the use of PSF-fitting methods and limiting us to aperture 
photometry. Aperture photometry was performed using the 
PHOT task in IRAF in all bands using the F435W+F550M 
detected catalogue. An aperture of 0.10" was used, and sky 
was sampled in a circular annulus of radius between 0.1 and 
0.15". Aperture cor rections were th en applied in accordance 
with those given in |Sirianni (2005) for the HRC and WFC 
bandpasses. While SBC/FI4OLP aperture correction s have 
been computed in the past, (e.g. iDieball et al.|[2005l ). they 
are not available in the published literature and our own 
aperture co rrections, computed with the TinyTim software 
l|Kristlll995l ). are included here in Table lAll Appendix [X] 
Since we have no a priori knowledge of the contirmum slope, 
we adopt the FI4OLP aperture correction for 13 = Q, using 
0.1" as in the other bandpasses. Since SB C/F122M contains 
the Lya line, no point-source photometry was performed in 
this bandpass itself. Instead, we perform aperture photome- 
try on the continuum subtracted Lya and Ha images in the 
same 0.1" apertures at the position of each cluster with no 
re-centering applied. This gives a direct measure of the Lya 
flux and equivalent width in the vicinity of each SSC. 

In order to estimate the properties of the SSCs, the 
same SED-fitting software described in Section 12.1.21 was 
applied to the aperture extracted fluxes. Age, E{B — V), 
photometric mass, etc. were computed, allowing us to com- 
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pare these properties with Lya in the immediate vicinity of 
each SSC. 



2.2 ESO - New Technology Telescope 

Haroll was observed during the nights of 18, 19, and 20 
September 2004, using the New Technology Telescope at 
ESO La Silla, as part of an observing run to obtain Ha, 
H/3, and [Olll] narrowband images for all southern targets 
in our EST Lya sample (Atek et al., in preparation). On 
the night beginning 18 Sept seeing was good, not exceeding 
1.2"for the duration, although thin cirrus prevents a direct 
calibration of these data. On the night of 19 Sept, obser- 
vational conditions were ideal: photometric and with sub- 
arcsecond seeing throughout. The final night was still pho- 
tometric but seeing deteriorated to > 2" and images were 
unusable for science purposes. The good-seeing data from 
the night of 18 was calibrated through secondary standard 
stars in the field, using data from the photometric nights 
of 19 and 20 Sept. Narrowband imaging was performed at 
Ha, H/3, and nearby continuum. Spectrophotometric stan- 
dard stars LDS749B, Feige 110, and GD50 were selected 
(on the grounds of their high spectral resolution: 1 — 2A, in 
order to resolve the stellar absorption features around Ha 
and H/3) from the catalog of JDke (1990), and were observed 
at regular intervals for the duration of each night in all fil- 
ters. Imaging observations were perfor med using both the 



ESO Multi-Mode Instrument (EM MI) JPekker et al 



and Super Seeing Imager 2 (SuSI2) (jP'Odorico et all 19981 ). 



19861) 



interchangeably. The instruments were used in 2 x 2 pixel 
binning mode to reduce the readout noise, providing a plate- 
scale of 0.1665" pix~^ and 0.130" pix~^, and a field-of-view 
(FoV) of 9.1 X 9.9" and 2.2 x 2.2", for EMMI-R and SuSI2, 
respectively. Table [2] summarises the imaging observations 
included in this article, lists the ESO filters used, and the 
total exposure times in each band. 

Data were first reduced by the standard routines in 
NOAO/IRAF: bias subtraction and fiat-field correction using 
well exposed sky- and dome-flats. Images were aligned us- 
ing the GEOMAP/GEOTRAN tasks and smoothed to the seeing 
of the worst seeing image using the GAUSS task. All HST 
bandpasses are aligned with the ATT frames and the PSF 
is degraded by convolution with a Gaussian kernel. Ha and 
H/3 were continuum subtracted to obtain the nebular emis- 
sion fluxes, accounting for contamination by [Nil] (Ha only) 
and underlying stellar absorption. [Nil] c ontamination is es- 
timated using [Nii]A6583A/Ha= 0.189 (|BergvaU fc OsthnI 
. 2002 ). Stellar absorption is estimated from the best-fitting 
Starburst99 spectrum at each pixel by modifying the spectral 
fitting code described in Sect. 12.1.21 The Starburst99 stel- 
lar libraries (|Martins et al.|[2005^ were built upon the latest 
model atmospheres and include full line-blanketing for all 
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stars and non-LTE effects for hot stars, thereby providing 
the best estimate of the stellar absorption features available. 
With the age of the stellar population we measure the equiv- 
alent width of H/3 using the same line and continuum wave- 
length windows as the Lick index l|Worthev et al ] |l994h . The 
windows used for Ha are of the same size as those for H/3, but 
scaled to Ha (i.e. the same windows x6563/4861 = 1.35). 
E{B — V) is generated from Ha/H/3, assuming a tempera- 
ture of 10,000 K and intrinsic line r atio of 2.86, us i ng th e 
extinction law of the SMC following [Gordon et aP l|2003h : 
iFitzpatrick fc Massal l|l988l ). 



circular region close to the source and avoiding other X-ray 
sources in the field. Source counts were grouped to have at 
least 20 counts per bin. SAS appropriate tasks were used 
to generate the distribution matrix and auxiliary response 
matrix files. 

On-axis observations of H aro 11 have also bee n obtained 
with the Chandra telescope ([Grimes et aLllioOTl l although 
their archive status is still proprietary. Quantities derived 
from that dataset are, however, compatible with those pre- 
sented here using off-axis archival observations (see Sec- 
tion [32J. 



2.3 X-Ray: Chandra and XMM-Newton 

Haro 11 happens to be located 14' away from the Cartwheel 
galax;y, of which X-ray observations have been obtained 
with both the Ch andra and XMM-Newton telescopes 
jWolter fc Trinchieril [2004). Respective total integration 
times in Chandra and XMM-Newton were 80 and 70 ksec. 
Regarding XMM-Newton, Haro 11 only falls within the field- 
of-view of the M0S2 detector, having unfortunately been 
missed by both the MOSl and PA'' chips. Both the Chan- 
dra and XMM-Newton datasets have been obtained from 
their respective archives. Since both telescopes operate with 
curved focal-plane configuration, the point spread function 
degrades severely with angular distance from the optical 
axis. For the XMM-Newton detection, the PSF is so dis- 
torted that all spatial information has been lost, but the 
object is bright enough that a one-dimensional spectrum 
can be extracted from the M0S2 data. In the Chandra ob- 
servation, Haro 11 lies on the SI chip: the object is clearly 
detected and we were able to perform a spatial analysis of 
the source (see Sect. 13. 2 |). 

The Chandra data were processed following the X-Ray 
Data Centre pipeline software. The level 1 events were re- 
processed using Ciao 3.4 task acis_process_events. The 
whole band (0.2-10 keV) image was smoothed applying the 
csmooth task. The image was smoothed with a minimal 
signal-to-noise of 3 using a circular Gaussian kernel. The 
source spectrum was extracted from a circular region of ra- 
dius 30". The background was obtained from the combina- 
tion of four circular regions located close to the source and 
avoiding other X-ray sources in the field. Source counts were 
grouped to have at least 20 cou nts per bin to allow mo di- 
fied minimisation technique ([Kendall fc Stuar3ll973h in 
the spectral analysis. Redistribution matrix and auxiliary 
response matrix files were generated. 

The XMM-Newton data were processed using the 

standard Sc ience Analysis System, SAS, v. 7. 0.0 

([Gabriel et al.[[2004t l. The most up-to-date files available as 
of January 2007 were used for the reduction process. The 
time intervals corresponding to high backgroun d events 
were removed using the method followed in Piconce lli et al.[ 
([2004h . The resulting exposure for the M0S2 data is only 
of 44 ks. No sign of pile-up was found in the M0S2 image, 
according to the epaplot SAS task. A visual inspection 
of the 0.2-10 keV image shows a distorted shape due to 
location of the source, close to the edge of the CCD. 
Therefore, no further analysis of the XMM-Newton image 
has been performed. The source spectrum was obtained, as 
for the Chandra observation, from a circular region with 
a radius of 30". The background was extracted from a 



3 RESULTS 

3.1 Lyman-alpha, UV and optical 

Figure [T] shows the results of the processed HST observa- 
tions. This panel of figures shows Haro 11 in FUV continuum 
and -B— band, continuum subtracted Lya and Ha, Wi^ya, 
and Lya/Ha. An inverted log-scale is used showing emission 
in yellow to black. Significant star-forming condensations 
2003) and the approximate kinematic centre 
19991 ) are marked in the _B— band image. Other 
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Ostlin et al. 



maps show the position of the SSCs and the Voronoi tessel- 
lation pattern. Plots at the bottom of the figure show the 
spatial variation of Lya, Wi^ya, and Lya/Ha along horizon- 
tal and vertical lines through the brightest Lya point, the 
positions of which are marked in the B— band image. 

We note that Ha emission loosely follows that of the 
FUV continuum (as would be expected from star-formation 
signatures) although there are a few differences. Notably, the 
brightest Ha knot (B) is not the brightest UV source (knot 
C). This could be an effect of age, dust, the density of the sur- 
rounding gas, or a combination of all three. The Lya images 
on the other hand, resemble neither the UV nor Ha. Knot C 
shows strong central Lya emission with two other patches of 
more diffuse emission to the north-west. Absorption is seen 
throughout the region of knot B and emission and absorp- 
tion can be seen varying on small (~ 30pc scales) around 
knot A, although overall the absorption dominates. A com- 
pact, isolated emission source to the south-east appears in 
the continuum subtracted image and is labeled knot D. More 
generally, however, the Lya emission appears in two distinct 
components: a central bright region and a surrounding dif- 
fuse, largely featureless, halo component. In the following 
analysis we identify these as the central and halo emission 
regions as two different physical components. 

Haro 11 has been observed through the low-resolution 
Short Wavelength Prime (SWP), 10x20" elliptical aperture 
of the lUE satellite. We obtained the lUE spectrum from 
the archive and have measured the lUE Lya flux to be 73.7- 
10"^* erg s~^ cm~^ with an equivalent width of 14. 8A. Pho- 
tometry of our Lya image in an aperture designed to match 
that of the /f/B reveals a flux of75.0T0~^^ ergs~^ cm~^ and 
equivalent width of 14. 3A. This flux is fully consistent with 
the photometric error on the lUE measurement of around 
10%^ 

[Kunth et al.l ([l998l ') presented Lya measurements ob- 
tained with the GHRS. This spectrum was obtained after a 
sequence of small angular maneuvers (SAMs), designed to 
find the peak in FUV surface brightness. Information about 
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Figure 1. Haroll as seen by HST. North is up and east is to the left. Figure sizes are 20 X 20", corresponding to a physical size of 
8.12 kpc square. Images are individually labeled. Main star forming condensations and the approximate kinematic centre are labeled in 
the S— band image and described in the text. The Upper right image shows the position of the star clusters - size is proportional to 
log(M) and colour shows age. The filtered Lya image has adaptively smoothed with ESO-MIDAS to give S/N = 20 for an additive 
noise model and overlayed contours are from FI4OLP. Log intensity scale is used throughout with high values showing up in yellow to 
black (i.e. black halos in the Wi^ya and Lya/Ho maps represent enhanced Lya). Lower plots show the spatial variation Lya, Wj^ya, and 
Lya/Ha in the rows and columns shown in the _B— band image. 
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Table 3. Lyo best-fit surface brightness profile parameters and 
integrated fiuxcs 



Component 


/LyQ.O 


[ " ] 


n 






Halo 


2.12 


1.05 


2.2 


76.4 


6.93 


Central excess 


47.8 


0.134 


1.0 


5.64 


0.512 


Total 








82.0 


7.44 


Measured total 








79.6 


7.22 



Notes: /Lya,0) in units of lO"^** erg s~^ cm^^ arcsec"^. 
^Lya in 10~^* erg s~^ cm~^. 
Li^ya in 10"*^ erg s"^. 

n is the Sersic index. Fluxes in the central and halo component 
fluxes are measured by integration of the Sersic function from 
zero to infinity and the Total flux is the sum of the two 
components. The Measured total is the numerical sum of all the 
flux measured in the the annuli. 

the SAMs is not available and the FUV morphology is so 
complex that it is not possible to reconstruct the maneuvers 
made by the telescope. In a 1.7 x 1.7" aperture centred on C, 
the continuum flux measured in the ACS frames is a factor 
of 4 brighter than that in GHRS, and apertures centred di- 
rectly upon A and B both show Lya in absorption. We are 
unfortunately not able to test our ACS' fluxes against this ob- 
servation, despite numerous efforts to predict the PEAK-UP 
behaviour: it's most likely the telescope centred somewhere 
between the three knots. 

The total emergent Lya line-flux measured in the image 
is 79.6 ■ 10~^* erg s~^ cm"'^, corresponding to a luminosity 
of 7.22 ■ 10*^ erg s~^ for our assumed cosmology. 

To eye the Lya image is so azimuthally symmetric that 
we deem it informative to present a radial surface bright- 
ness profile of continuum-subtracted Lya, shown in Figure (2] 
This was generated by integrating the flux in concentric an- 
nuli of width 2 pixels, centred on the brightest central pixel 
of knot C. The central and halo components identified in 
the LyQ map are clearly visible in the surface brightness 
profile, with the central component taking over at around 
0.5". Since they are so clearly independent (visually resem- 
bling a bulge-|-disc morphology) we fit Sersic functions to 
the halo and central excess (i.e. total-halo) components us- 
ing the formula 

/A,R = /A,oexp[-(i?/i?o)'^"] (1) 

where /a,o is the central surface brightness, Ro is the effec- 
tive scale length, n the Sersic index. Integration from zero 
to inflnity then gives back the total flux in each component. 
The results of the proflle fltting in the two components can 
be seen in Table O The regions chosen for fltting in the two 
components were 0-0.5" and 1-7" for the central region and 
halo, respectively. 

We have also Voronoi binned the individual images. 
Since Voronoi binning conserves surface brightness across 
each spaxel, we can directly compare Lya surface brightness 
with other fluxes or properties without a loss of reliability 
in lower surface brightness regions. With binned images in 
hand, we present the Lya surface brightness and Lya/Ha 
as a function of various other parameters in Figure [3] Each 
point represents the surface brightness in each resolution 
element, although the elements themselves vary from single 




0.1 1 10 

Radius [ arcsec ] 



Figure 2. LyQ surface brightness profile. The vertical dashed 
line indicates the radius at which the annuli intersect with the 
nearest edge of the SBC detector. Green and red lines show the 
best-fitting Sersic profiles to the halo component, and the central 
excess, respectively. The blue line shows the sum of the surface 
brightnesses in the two components. 



pixels to 1600-pixel (1 □") spaxels. Also, since fluxes are pre- 
sented on log-space, only resolution elements that are posi- 
tive in Lya are shown, eliminating regions of absorption. It 
is not possible to estimate global or integrated values from 
these scatter-plots; histograms of the distributions are pre- 
sented above and to the side of each scatter plot, in order to 
better give a feeling for this. The left-most two plots show 
how the Lya emission and line ratio correlate with Ha and 
the age of the young stellar population. The Lya vs. Ha plot 
is again suggestive of two distinct pixel components. One at 
lower Ha surface brightness (< lO^^erg s~^ cm~^ arcsec"^) 
in which Lya appears largely independent of Ha. This fea- 
ture extends over approximately three orders of magnitude 
in Ha and the Lya/Ha plot shows the emission is always 
above the value that would be expected from recombination 
theory. Above the quoted Ha surface brightness, the distri- 
bution changes and an upper limit appears to be set by the 
recombination value. The line ratio plot clearly shows that 
Lya photons are emitted from regions where the Ha flux 
would suggest that little nebular gas is present. A similar 
effect is also seen in the second plot where Lya is compared 
with the age of the stellar population and there is evidence 
for substantial amount of pixels showing Lya emission with 
ages in excess of 100 Myr; much too old to produce nebular 
emission. The right-most two plots show how Lya correlates 
with the E{B — V) as determined from fitting the stellar 
UV continuum {E{B — V)*), and continuum slope itself (/3). 
Since /3 in part refiects E{B — V), these two plots correlate 
somewhat in shape, although /3 is also dependent upon age. 
The discrete values of E{B — V)t from the fltting grid are 
apparent. Again the distribution signifles a non-continuous 
distribution of pixels: in both cases Lya can be seen in emis- 
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Figure 3. Resolution element scatter plots of Lya surface brightness [upper] and the Lyo/H« line ratio [lower] against various measured 
values. Far left: Ha surface brightness. Centre left: Age. Centre right: Dust reddening E[B — V). Far right: FUV continuum slope /3 
between 1500 and 2200A. Red lines show the recombination ratio of Lya/Ha in the case B approximation and the Lyo=Ha line. Green 
lines in the E[B — V) plots show how dust attenuates Lyo for various intrinsic fluxes [upper) and how it reduces the line ratio [lower). 
The resolution elements themselves vary in size and are not representative of absolute fluxes. The histograms above and to the side of 
the scatter plots show the distribution of the values, in order to give some impression of how the various parameters contribute globally. 



sion from regions with E[B — V) extending up to 1, and 
(3 showing a very red slope. In these regions of the plots, 
emission is again found almost exclusively above the recom- 
bination value. At lower E[B ~ V) [< 0.3) and /? (< -1) 
the pixel distribution changes and broadens significantly in 
both Lya and Lya/Ha. The Lya surface brightness tends 
to increase in these regions although the bulk of the pixels 
show line ratios significantly below the recombination value. 
While brighter overall, these emission regions show strong 
attenuation of Lya. 

Figure U shows how FLya correlates with gas- phase in- 
terstellar extinction, E{B — V)is, as determined from the 
Balmer decrement and is partly similar to the right-most 
plots in Figure O This provides a much less biased esti- 
mate of dust extinction compared to the approximation ob- 
tained by spectral modeling of the stellar continuum. Ta- 
ble|3]shows the values of E[B — V)is obtained by measuring 
Ha and H/3 in 1" circular apertures centred on knots A, B, 
an d C, together with the values obtained spectroscopically 
bv lVader et al.l l|l993l ). With the spatial sampling in HST 
images degraded by two orders of magnitude to match the 
ground-based seeing, we cannot resolve the fine details of 
the Lya emission and the surface brightness range is signif- 
icantly reduced. -pLya is here presented in linear space and 
now shows regions of absorption, most notably around knot 
A. This distribution shows a structure partially resembling 
that of the Lya vs. E{B — 1/)* plot: a component of difi'use, 
low surface brightness Lya emission is seen extending be- 
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Figure 4. FLya vs. interstellar E[B - V) [E[B - V)is) as de- 
termined by the Balmer decrement. Regions of knots A, B, C, 
and the diffuse emission component are represented in different 
colours and marked in the figure. 



yond E[B — V)is^ 1, with a brighter Lya emission region 
(knot C) where E[B — V)is shows a large spread, with a 
mean value of 0.48. 

In total we find 140 star-clusters. The spatial distribu- 
tion, together with the ages and masses of these clusters can 
be seen in the upper-right panel of Figure [1] The age and 
mass distribution of the clusters can be seen in Figure [5] The 
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Table 4. E{B-V)is as measured in the NTTHa and H/3 frames, 
averaged in circular apertures centred on the 3 main star forming 
condensations. Values obtained by Vader et al. (1993) are shown 
for comparison. 
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Figure 5. Upper. Age distribution of the star clusters; Lower. 
Mass distribution. Knot C is dominated by a single SSC with a 
mass of 2.2 X 10® Mq, with an age of around 2 Myr. 



decrease in the mass distribution at < 10'' M© is the result 
of incompleteness in the less luminous clusters. While there 
is much of interest in these data, most of it is outside the 
scope of this article and will be presented in a forthcoming 
publication (Adamo et al., 2007 in prep); here we simply 
deal with the relation between a few basic properties and 
Lya. 



3.2 X-ray: XMM and Chandra 

In the upper panel of Figure |S] we show the Chandra image 
of the full energy range (0.3-8 keV) adaptively smoothed 
(see Sect. 12. 3|) . Continuum contours from the ACS/FI4OLP 
observation have been overlaid to show the location of the 
different knots. The lower panel in Figure |6] shows the X- 
ray surface brightness profile, in the 0.2-2 keV range com- 
pared with the corresponding profile of the off-axis PSF of 
ACIS-Sl at 1.5 keV. The plot shows hints of extended emis- 
sion visible in the 20-40 px region, i.e. ~10-20 " (1 pix = 
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Figure 6. Upper. Haroll as seen by Chandra over the complete 
0.3-8keV spectral range. Image is 40" x25". Overlaid are the FUV 
contours as seen by HST/FI4OLP. Lower. Chandra X-ray source 
surface brightness profile over the 0.2-2 keV range together with 
the corresponding PSF at 1.5 keV. Datapoints consistently exceed 
the PSF between ~ 10 and 20", showing the source to be resolved 
and spatially extended. 



0.492"). The Chandra image shows a central unresolved core 
region and a weak, diffuse hot gas component, emitting over 
an area ^ 30" in diameter. The diffuse component appears 
more extended in the north-south direction than in east- 
west. However, no firm conclusion can be reached as the im- 
age is largely distorted due to the ~ 14' offset from the Chan- 
dra nominal on-axis location. In fact, the Chandra on-axis 
images of Haroll reported by G rimes et al.. (.2007. ). shows 
different emission knots following the optical/UV emis sion 
structure. The brightest peak of the lGrimes et all (|2007l ) im- 
age coincides with knot B, while the maximum emission on 
our image is located closer to knot C and more than 5" from 
knot B. A plausible explanation for the observed mismatch 
could be variability, although the spectral analysis, fluxes, 
and luminosities of both observations are in good agreement 
(see below). Therefore, the more likely explanation could be 
the intrinsic Chandra astrometry uncertainty, although we 
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Figure 7. Chandra spectral extraction in the 0.2-10 keV range, 
the best fit model, and the residuals. Model parameters can be 
found in Table [5] 

note that this is only of the order of 2" for off-axis observa- 
tions. 

Spectral analysis was performed on both Chandra and 
XMM-Newton data, in order to derive the spectral shape of 
the source, the X-ray luminosities and, local photoelectric 
absorption column density of Hi. The observed spectrum 
was corrected for Galactic absorption, w hi= 2.39- 10^" cm^^ 
iKalberla et al.ll2005l : lBaiaia et aLlbOOSl ). The Haro 11 spec- 
tral emission is complex: both power law and thermal single- 
component emission models provide unsatisfactory fits to 
the data, with values equal to 2.2 and 3.5, respectively. 
The best-fit model consists of a comb ination of a power-law 
and mekal model (|Mewe et al.l 11986^. Statistically indistin- 
guishable results were obtained applying a double thermal 
mekal model. However, the temperature of the hottest ther- 
mal model is 9^2 keV. This high value is difficult to explain 
in a standard physical scenario. We have also derived an up- 
per limit for the absorption equivalent Hydrogen column of 
6 X 10^° cm~'^. Figure [7] shows the X-ray spectrum together 
with the best fit model and the corresponding residuals. The 
values of the parameters and the goodness of the fit are 
listed in Table [S] The corresponding total X-ray fluxes and 
unabsorbed luminosities are shown in Table [6] Even though 
the effective area of XMM-Newton is larger than that of the 
Chandra, the M0S2 spectrum of Haro 11 resulted in lower 
signal-to-noise. This is mainly due to the short real expo- 
sure time (44 ks) after the high background filtering, com- 
bined with the fact that the source is only visible in M0S2, 
for which the effective area is 3 times lower than that of PN. 
However, the source MOS spectra was analysed, confirm- 
ing the results obtained with Chandra. A single-component 
model (power-law or thermal emission) failed to reproduced 
the data. The best fit model consists also of a power-law 
with r ~ 2 and a mekal component with kT ~ 0.4. All 
the parameters derived from the XMM-Newton spectrum 
are compatible within the error with those from the Chan- 
dra analysis. We have measured a somewhat lower value of 
the high energy fiux, of the order of 65±20%. The quanti- 
ties derived from our spectral analysis are compatible with 
th ose derived fr o m on -axis Chandra observations presented 
in lGrimes et all (|2007l '). 
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Table 5. Best-fitting parameters of X-ray spectral model. 



Parameter 


Value 


Unit 


Intrinsic nui 


< 6 X 10^0 


cm~-^ 


Power-law photon index 


1.72 ±0.18 




Mekal kT 


'-'■0^-0.04 


keV 



X^: 95.30 using 88 PHA bins. 

Reduced x^: 1-08 for 88 de erees of freedom. 

Table 6. Observed X-ray fluxes and unabsorbed luminosities. 



Band [ keV ] 


Flux 


Luminosity 


0.5 to 2 
2 to 10 


q 7+0.4 
^- ' -1.3 

13.3«;^ 


1 20+"-°* 


Total 


23.0t^;* 


2 2fi+'''l'' 



Flux in 10 erg s ^ cm ^ 
Luminosity in lO**^ erg s^^ 

4 ANALYSIS AND DISCUSSION 
4.1 Lyman-alpha morphology 

Figure [T] shows some of the imaging results. Visual inspec- 
tion clearly reveals that Lya morphology is not represen- 
tative of either Ha or the FUV continuum. The overall 
Lya morphology exhibits a central bright emission source 
(centred on knot C), with a diffuse, largely symmetric halo- 
like component with lower surface brightness. The emission 
around C is more extended in the north-south directions 
compared to east-west, and appears as two fan- like struc- 
tures of size ~ 1". This structure is not seen in E{B — V) 
maps and may be due to a filamentary Hi structure. Typi- 
cally Why a and Lyo/Ha are low in the central regions, where 
continuum and H« are strong, but increase dramatically in 
the outer regions. Some absorption spots are seen in the 
regions around knots A and B. The most likely interpreta- 
tion of this overall morphology is that photons are escaping 
directly from the central regions, either through winds or 
ionised holes in the ISM, while the halo component arises 
from resonance scattering of Lya. 

Knot P, the small blob o f Lya emission found in the 
analysis of (|Kunth et al.ll2003l ') is clearly visible towards the 
south-west. Knot T> is not visible in any of the continuum 
bands or the Ha on-line observation and it appears to be a 
real source of Lya photons devoid of continuum. 

The radial surface brightness profile presented in Fig- 
ure [5] shows the fiux in the halo component decaying ap- 
proximately exponentially with radius, but peaking sharply 
in the centre. The fan-like structures described above are 
contained within the central surface brightness peak. Fitting 
Sersic functions to these two independent components and 
integrating under the best-fitting formulae allows us to mea- 
sure the relative contribution of the components of central 
and halo emission, and shows that while surface brightnesses 
are sigrnficantly lower, ~ 90% of the photons are emitted via 
the photon diffusion mechanism. This is comparable to the 
value of ~ 70% in the halo com ponent of ESO 338-IG04 pre- 
sented in iHaves et al.1 (120051 ). although for that study the 
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estimate was made by masking the images instead of profile 
fitting. The two galaxies are comparable in FUV luminosity, 
both low-metallicity, and both exhibit complex dynamical 
and morphological structures, and the Lya output of both 
are clearly dominated by the halo mechanism. 

4.2 Lyman-alpha production, escape and 
regulation 

In the E[q image we measure a total line flux of 333 x 
10"^** erg s~^ cm~^, identical to that given in lSchmitt et al.l 
l|2006l ). which is not corrected for Galactic extinction. Galac- 
tic extinction correction (jCardelli et al.l Il989l ) gives a in- 
trinsi c Hq flux o f 342 x 10"^* erg s~^ cm~^, consistent 
with lOstlin et al.l (|l999l ). The total Lya production esti- 
mated from the case B recombination gives a value of 
2975 X 10~^* erg s""*^ cm~^. Adopting the Lya flux we 
obtained by integrating the (Galactic-extinction corrected) 
surface brightness profile we obtain a Lya escape fraction of 
just 2.7%. This calculation neglects the fact that Lya can 
be enhanced relative to Ha by shock excitation. Were there 
a significant contribution to Lya as a result of a shocked 
ISM, this escape fraction will be further reduced. 

Plots of the Lya surface brightness in each binned el- 
ement, compared with various other quantities, are seen 
in Figure [3] Most importantly, these plots show how the 
Lya surface brightness spatially compares with Ha and 
the dust distribution. The left-most plot shows Lya vs. 
Ha; the red lines indicate Lya=Ha line and Lya/Ha=8.7 
(case B recombination). Comparing Lya to Ha also identi- 
fies two independent components of Lya emission. Firstly, 
at low Lya and Ha surface brightness (< 10"^'' and 
10~^^ erg s~^ cm~^ arcsec"^ for Lya and Ha, respectively), 
the Lya to Ha ratio seems to always exceed that predicted 
by recombination (see Lya/Ha spatial distribution in Fig- 
ure [T]). Again, shock excitation could enhance the produc- 
tion of Lya relative to Ha, and it's plausible that some of 
the Lya is produced in a shocked ISM. Indeed, this is an 
oft-proposed mechanism for Lya production in high-z Lya- 
blobs (LABs), resulting in large-scale Lya emission with lit- 
tle or no associated continuum. However, without observa- 
tions in other emission lines the exact mechanisms at play 
in LABs are very difficult to test, although tests may be 
performed for targets at low-z observed with HST. Opti- 
cal sp ectros copy of Haro f f perf ormed bv iHeisler fc Vadeij 
l|l994 ) and iKewlev et al.1 (120011 ) consistently find spectra 
consistent with pure photoionisation nebulae. On the other 
hand, these spectra were extracted in large apertures centred 
on the knots, avoiding difi'use Lya emission regions, and the 
answer is not fully conclusive. It is possible that future high 
resolution spectroscopy (e.g. with a repaired STIS) may be 
able to test this directly. Given the resonant nature of Lya 
a more likely interpretation of this Lya-Ha decoupling is 
that it arises from resonant scattering of Lya photons. The 
Lya/Ha histogram in Figure |3]confirms that over 50% of the 
pixels are in the component with super-case B line ratios. 
The fact that these points are here most likely indicates that 
a significant number of Lya photons must have travelled far 
from their point of origin, undergoing numerous resonance 
scatterings. 

At higher surface brightnesses (> 
10~^^ erg s~^ cm~^ arcsec"^ for Ha), the distribution 



clearly changes. Here the recombination value seems to 
set an upper limit on the Lya/Ha ratio and Lya is never 
stronger than expected. The top edge of this distribution 
must result from Lya photons escaping the galaxy directly 
either due to ionised holes in the ISM cleared by the UV 
radiation field or through a rapidly outflowing neutral 
medium. The distribution fans out in this component 
and a large spread in Lya surface brightness is seen which 
results from either resonance scattering selectively removing 
Lya photons from the line-of-sight (from where they may 
be re-scattered and emitted in the halo component) or 
destruction by dust. The central emission region also shows 
absorption structures to the east and west of knot C (see 
middle-right panel of Figure [1} , where this effect may be 
strongest. 

The resonant decoupling of Lya photons is also re- 
vealed in Figures |3] and U where Lya is compared with 
E{B — V) from fitting to the stellar continuum {E(B — 
V)*), Ha/H/? {E{B - V)is), and /3. If Lya were a non- 
resonant line a declining relationship would be expected 
between Lya and E{B — V). The bulk of the Lya out- 
put is in the halo component (Lya surface brightness < 
10^^'' erg s~^ cm~^ arcsec"'^; Figures [2] and |3} ■ Below this 
level, there is essentially no trend of Lya with E{B — V)t. 
Most importantly, when Lya is seen in emission from re- 
gions with E{B — V)t> 0.2 it is almost entirely with line ra- 
tios that exceed the recombination value. This uncorrelated 
faint emission component is also seen in the E{B — V)is plot 
where the distribution of diffuse -pLya pixels (shown in cyan) 
is constant over the whole range of E{B — V)is. Despite the 
low dust content of knot A (mean E{B — V)is~ 0.2), Lya 
is seen only in absorption in this region. In contrast, the 
brightest Lya region (knot C) shows Lya only in emission 
with measured values of E{B — V)is as high as 1, although 
the regions A, B &i C have sizes comparable to the seeing 
disc in the ground-based images and therefore the scatter at 
each knot may be an artefact (the mean value of E{B — V)is 
in knot C is 0.48). Lya reaches its brightest flux levels where 
E{B — V)i, is small which is symptomatic of the contrasting 
escape mechanisms, but it must be noted that this only ac- 
counts for around 10% of the overall Lya output. That we 
see such strong absorption from knot A but emission from 
knot C seems to indicate that the static Hi column density in 
the line-of-sight to knot A must be significantly greater than 
to knot C. Knot A does contain several massive clusters, the 
ages of w hich are consistent w ith the high mechanical lu- 
minosity (iLeitherer et al.lll999l ) that would be required to 
drive out the neutral ISM, although it seems that either 
the Hi covering is too great, or the clusters are too young to 
have accelerated the surrounding medium. Unfortunately we 
have no direct information on the Hi column density or the 
mechanical energy return that would be required to drive a 
superwind. 

That we see Lya from resolved regions with E{B — V)> 
0.5 is a striking result. For a simple dust slab this would 
correspond to the loss of more than 99.9% of Lya photons, 
before scattering is considered. Since in the halo regions we 
typically see Lya/Ha above the recombination value we can 
infer that radiative transfer effects are at play in the prop- 
agation of Lya photons, and scattering events have caused 
these photons to travel significantly increased path-lengths. 
Again, knowledge of the geometry and kinematics of the neu- 
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tral and ionised media would be required to infer the num- 
ber of scatterings, mean free paths, and therefore the actual 
path-lengths involved. We can, however, identify two pos- 
sible extremes of the photon transport mechanism. Firstly, 
they may have diffused through largely static Hi gas with 
transport dominated by multiple rare scatterings into the 
wings of the distribution function. In this case, the neutral 
medium must be pristine, dust-free, and have never mixed 
with star-forming regions. Secondly, and much more likely, 
much of the propagation may have been through paths in 
the low-density ionised ISM with scattering/reflection events 
occurring at the surfaces of dense neutra l clou ds. In the mul- 
tiphas e ISM models of iNeufeldl (|l99ll ) and [Hansen fc Olil 
120061 ) ■ such a geometry shields Lya from the dust that is 
embedded in the Hi, with reflections preserving Lya by pref- 
erential propagation in the ionised ISM. Transport effects in 
these models suggest that equivalent widths may be signif- 
icantly boosted, even globally. The global Lya equivalent 
width and line-ratio are certainly not boosted, and what we 
are seeing may be preferential Lya destruction in some re- 
gions and relatively unhindered transport in others. Where 
the bulk of the Lya destruction is taking place, we cannot 
say. Either way, the transport effects result in the super- 
position of Lya photons on top of dusty and/or old stellar 
populations. 

A two dimensional kinem atic study of Haro 11 has been 
performed in the Ha line by lOsthn et aP (| 19991 . 1200 ih us- 
ing the C/G^iS Fabry-Perot (FP) interferometer on ESO's 
3.6m telescope. These data demonstrated the kinematic cen- 
tre to lie approximately midway between knots A and C as 
labeled in Figure [T] with the major axis lying in the NW- 
SE direction with the NE side toward knot B receding. In 
the central 2"(lkpc) iso- velocity contours are shown to be 
densely packed and show a steep rise in rotational velocity 
in the direction of knot B, indicating a non-relaxed system 
or possible second kinematic component. While we expect 
even the diffuse Lya emission to be in part regulated by the 
kinematic structure of the ISM, we are not able to see any 
correlation between Lya and velocity shifts or shear in the 
ISM. The resolution of the CIGALE data is not sufficient 
for us to compare small-scale Lya emission/absorption fea- 
tures with the kinematic structure of the ISM, although on 
larger (~kpc) scales no correlation is apparent. Moreover, 
turbulence or disruptions in the ionised medium does not 
necessarily imply the neutral ISM has been accelerated. 

In summary we see clear evidence that, while Lya pho- 
tons must be produced in the same regions as Ha, we do not 
observe them to be emitted from their production sites, and 
the bulk of the Lya photons travel significantly enhanced 
path lengths. We see no relationship that would indicate 
Lya correlates spatially with the dust distribution. Unfortu- 
nately we do not have information about the Hi distribution 
or its kinematic structure. 



4.3 Lyman-alpha and the SSC population 

The age distribution of the super star cluster population are 
shown in Figure [S] Figure [8] shows the equivalent width of 
Lya in the immediate 0.1" of the super star clusters, com- 
pared with the respective age of each cluster. Lya equiv- 
alent widths are measured from the continuum-subtracted 
and continuum-only image in fixed apertures. 




Age [ Myr ] 

Figure 8. Distribution of Lyo equivalent width in the immediate 
vicinity of the each super star cluster compared with the cluster 
age. 



iTenorio-Tagle et all l| 19991 ) and lMas-Hesse et all (|2003l ) 

have presented evolutionary studies of Lya emission from 
starbursts which qualitatively predict the evolution of Whya 
with age. The Lya sequence is more complex than that of Ha 
due to the resonant trapping of Lya photons in Hi resulting 
in a need for mechanical feedback and/or a higher ionisation 
fraction to create a medium optically thin to Lya. The model 
is such that at very young ages (< 1.5 Myr) the Lya feature 
is expected in absorption due to dense, static Hi medium. 
After some time, the cluster is able to ionise this medium 
and drive out an expanding shell resulting in a pure emission 
feature between 2.5 and 5 Myr. As the expanding shell re- 
combines, the Lya emission feature develops a P Cygni-like 
profile which, with continuing recombination removes suc- 
cessively more of the emission until, at around 8 Myr, the 
Lya feature is again in absorption. With such a large num- 
ber of SSCs (effectively mini galaxies) and Lya maps, this 
evolutionary sequence becomes something we could possibly 
examine within our dataset. However, there are a number 
of factors which could mask such a relationship, particu- 
larly in a dataset consisting of imaging observations only. 
Firstly, when the feature is P Cygni some or all of the emit- 
ted fiux is cancelled by the absorption profile, resulting in 
termination of the observed Lya emission phase at earlier 
times than predicted. Secondly, geometry also enters into the 
model and, due to the ionisation structure being cone-like, 
the orientation must be favourable in order to observe the 
Lya emission phase. In our dataset, this will result in a de- 
crease of Whya or complete absence of Lya from some Lya 
emitting clusters. Thirdly, due to the diffuse emission, Lya 
can be seen from essentially any region of the galaxy. How- 
ever, due to the high continuum levels around SSCs, strong 
contamination is not expected and will also be accounted for 
by background subtraction. Finally, a general spread can be 
expected due to variations in the initial cluster mass. Hi dis- 
tribution, mechanical feedback requirements, and internal 
absorption that is not local to the cluster itself. What we 
could possibly expect is negative Whya at the earliest times, 
becoming positive and peaking between 2.5 and 4 Myr, and 
then declining again to the background levels at around 6 to 
8 Myr. Spread is expected toward lower (and negative) Whya 
due to orientation effects and the Hi covering. Moreover, the 
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Table 7. Integrated star-formation rates in Haroll. The VLA 
data have been taken from Condon et al. (1998). 



Tracer 


SFR 


Source 


Calibration reference 




1.4 GHz 


5.6 


VLA 


Condon (1992) 




FIR 


18 


IRAS 


Rosa-Gonzalez etaL 


(2002) 


Ha 


24 


HST 


Kennicutt (1998) 




1527A 


19 


HST 


Rosa-Gonzalez et al. 


(2002) 


0.5-2 kcV 


22.oj;«-« 


Chandra 


Ranalli et al. (2003) 




2-10 keV 




Chandra 


Ranalli et al. ("2003) 




SFR [ Mq yr 











noisy nature of the continuum-subtracted Lya map further 
complicates very local equivalent width measurements. We 
consider only the sub-population of clusters with age below 
40 Myr. Figure [8] does show a distribution of peak Wi^ya vs. 
age that is quantitatively consistent with the evolutionary 
models. For the lowest ages Lya is seen only in absorption 
or around the zero-level. For ages between 2 and 4 Myr some 
SSCs show local Lya emission that could possibly represent 
the emission phase, although there is some spread. At late 
times (> 8 Myr) Lya is only seen in absorption or very close 
to zero. 

4.4 X-ray information and star formation rates 

The Chandra image (Figure |6]) shows that most of the X- 
ray photons have been emitted in a region surrounding knot 
C. The radial profile of X-ray emission is largely feature- 
less. Comparing the radial profile with the Chandra off-axis 
PSF (Figure [S| yields some hints of extended emission up to 
around 10-20" . The X-ray emission arises from the same re- 
gion around knot C where Lya is brightest, suggesting that 
Lya emission from this region might be the result of a per- 
turbed and outflowing ISM, accelerated by the release of me- 
chanical energy from the SSCs located at knot C. However, 
it is worth noting that the analysis of an on-axis Chandra 
observation of Haroll (.Grimes et al.ii2007i ) indicates that 
the emission is resolved in several knots, the brightest one 
coinciding with knot B. 

We have measured an upper limit for the equivalent 
Hydrogen absorption column of 6 x 10^" cm~^. We want 
to stress that consistent fits could also be found with lower 
columns: the F-test gives only a probability around 90% for 
the fits using this value of the column density, so that the 
present data do not allow us to constrain the amount of neu- 
tral Hydrogen in front of knot C . In any case, column den- 
sities around 2 x 10^" cm~'^ were obtained bv iKunth et al.l 
l|l998t ) from Voigt profile fitting of the Lya absorption 
pro file in the GHRS spectrum (hhi = 2.5 ■ 10'^'' cm~^), 
and IVerhamme et aL 

I ([2006) note that Voigt profile fitting 
may systematically underestimate nm- It is interesting to 
note that Haro 1 1 showed the largest Hi column density 
in the 4 galaxies of iKunth et all l|l998h with Lya emission 
(log(nHi) ~ 19.7 — 20.4 cm"'^). The derived temperature 
of the hot diffuse component is ~ 8 ■ 10''K, a typical value 
also found for other local dwa rf-starbursts like He 2-10 and 
NGC4449 (|Grimes et al.ll2005l ). 

In Table [7] we present star-formation rates (SFR) de- 
rived from X-rays, other observations in this study, and pre- 



vious observations from the literature. Within the errors, 
star formation rates derived from FUV continuum. Far In- 
frared (IRAS FIR, corrected to the 1-1000 \im range accord- 
ing to iHelou et al.l (.1988)), Ha and X-ray luminosities are 
consistent with a SFR ~ 20 - 25 M© yr"\ Such a SFR is 
rather high and could explain the presence of high-velocity 
winds and disrupted interstellar medium, driving the Lya 
escape and pushing hot and neutral gas to a large radial 
exte nt. Compared to the s ample of Lya emitters at 2 ~ 3 
from lGronwall et all (|2007l ). Haro 11 appears at the top end 
of the SFR distribution they derive from the FUV contin- 
uum (within the range 1-40, with an average value around 
5 M0 yr~^). This implies that Haroll is a very actively 
star-forming galaxy in the local Universe, comparable to the 
objects being analysed at high redshifts and the results ob- 
tained from its analysis may therefore be directly applicable 
to more primeval systems. 

The fact that the SFRs derived from the X-ray emis- 
sion are pretty consistent with the values derived from FUV 
continuum. Ha and FIR emission implies that the X-ray 
emission should be driven by the powerful star formation 
processes taking place in Haro 11, with no hint on the pres- 
ence of a low luminosity active nucleus at the core of knot 
C. Indeed, the Lo.5-2keV /Lfir ratio of 2.6 ■ lO"'^ i s very 
close to the average value found by iRanalli et al.l (|2003f ) 
for a sample of star-forming galaxies at redshifts in the 
ra nge 0.2-1.3. A preliminary analy sis based on the models 
of lOti fc Mas-Hessd (|2007. in prepf) shows that this ratio is 
consistent with a powerful star formation process dominated 
by SSCs with age below 5 Myr, with a standard efficiency 
of around 10% in the transformation of released mechanical 
energy into thermal X-ray emission. 

On the other hand, the radio luminosity yields a signifi- 
cantly lower value of the star formation rate when compared 
to the other tracers. We believe this is due to the youth 
of the star formation episodes in Haroll. The calibration 
used assumes a steady star-formation process with radio 
contributions from both thermal and non-thermal compo- 
nents (synchrotron emission from supernova remnants). Fig- 
ure [5] shows that most of the SSCs are younger than around 
3.5 Myr, so that the number of supernovae expected from 
the slightly older SSCs is significantly smaller. The low radio 
emission could therefore be explained by this deficit in the 
non-thermal component. 

The relative hardness of the X-ray continuum (photon 
index ~ 1-7), indicates that the hard band, contributing 
~ 50% of the X-ray output, might be dominated by high- 
mass X-ray binaries (HMXB). Since HMXBs become active 
after 3-4 Myr of evolution, we believe that these sources 
could be already present in the SSCs with age between 4 
and 10 Myr. Indeed, an SSC of mass ~ 1.1 x 10*^ Mq with 
an age of 8.2 Myr identified near knot C could be already 
hosting a significant number of HMXB. The observational 
properties are therefore consistent with the distribution of 
ages shown in Figure [5] old enough to have produced the 
number of HMXBs seen in the hard X-ray bands, but still 
young enough not to have generated so many SNe that the 
radio luminosity would be dominated by synchrotron emis- 
sion. 
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4.5 Lyman continuum production and emission 



FUSE observations (jBergvall et al.l l2006h have shown 
Haro 11 to be an emitter of LyC radiation, with an observed 
flux density at 900A (/9oo,e) in the 30"x30" aperture of 
1.1 X 10~ ^^ erg s~^ cm~^ A ~^. Subsequent re-analysis of 
the data ijGrimes et alJ [20071 ) place /goo.e at a lower value 
of 0.23 X 10~^* er g s~^ cm~^ altho ugh the answer is 
not yet conclusive. iBergvall et al.l (|2006l ) estimated the in- 
trinsic LyC production by summing the measured LyC leak- 
age and an estimate of the ionising flux computed from Hq 
and recombination theory. This translated to an escape frac- 
tion of between 4 and 10% assuming a Salpeter IMF, al- 
th ough the higher valu e is reduced to ~ 2% in the analysis 
of[G rimes et al.l(|2007l ). However, this method of estimating 
/goo.o is difflcult to calibrate for a given wavelength since Hi 
is ionised by all photons with A < 912A, not just those at 



As described in Section [2TL2] our software also outputs 
maps of the intrinsic /goo, based upon the normalisation of 
best-fitting population synthesis models. This relies upon 
stellar continuum and allows us to predict directly the value 
of /goo.o at each pixel, thereby giving a different, and possi- 
bly more robust handle on the ionising photon budget than 
that estimates based upon Ha alone: Ha estimates provide 
and upper limit on the stellar ionisation budget since shock 
excitation of Ha is neglected, while stellar continuum fit- 
ting cannot account for ionising photons that arise from 
non-thermal processes. However, in light of the tight cor- 
relation between between X-ray determined SFR and SFRs 
determined by other methods, we do not expect a significant 
non-thermal contribution. 

Figure [9] shows a map of the ionising photon produc- 
tion, zoomed in on the central star-forming regions (i.e. the 
sites of LyC production). As expected, the LyC morphology 
itself resembles that of Ha, with knots A, B, and C show- 
ing up as strong production sites. Knot B is the brightest of 
the three in LyC, consistent with indications from the Ha 
that it is the most intensely ionising source. An apparent 
dust lane is seen north-west of B which is seen in the FUV 
continuum and optical bands, and may represent the fail- 
ure of the implementation of 'dust screen' extinction laws to 
recover all the intrinsic 900A flux when dust is distributed 
along the line-of-sight. Since the method of computing /goo 
relies on the normalisation of model spectra, /goo can never 
be negative and pixel-to-pixel noise can introduce a spu- 
rious positive contribution to the integrated flux. Here we 
deal only with the Voronoi-binned images in order to better 
control the background noise. 

bumming the 900A flux we obtain a total value of 
/goo.o = 12.3 X 10~^^ erg s~^ cm~^ A~^. This translates into 
a 900A escape fraction, /esc, of 9 %. This is towards th e up- 
per end of the range presented in IBergvall et al.] (|2006| ) but, 
again, neglects 900A photons of non-stellar o rigin. The up- 
per lim it on the escape fraction published bv lGrimes et al.l 
120071 ) ■ whic h compared only wit h the higher estimate of 
/goo.o from (jBergvall et all (2006l l. remains unchanged at 
around 2%. 

FUSE data contains no spatial information about the 
LyC leakage. Both LyC and Lya are strongly affected by Hi 
covering but the nature of the effect is different for the re- 
spective photon frequencies: LyC photons are absorbed but 




Figure 9. Production map of the flux density at 900A. North 
is up and east to the left. Image dimensions are 12 X 12", and 
concentrates only on the central star-forming regions where LyC 
photons can be expected to be produced. 



Lya photons are resonantly scattered, resulting in spatial re- 
distribution. Since LyC photons do not resonantly scatter, 
their emission most likely traces their production sites. The 
spatial redistribution of Lya can, to some extent, be mapped 
by our imaging technique. Only the central Lya emission 
region shows Lya/Ha around case B level and hence from 
knot C, Lya is likely to be leaving the galaxy directly due 
to an outflowing neutral medium, or cleaned ionised holes. 
The fact that Lya may directly escape the galaxy from knot 
C is suggestive that some LyC photons may also be doing 
so and, in the limit that all the LyC radiation leaks from C 
we obtain fasc = 1/3 for LyC at knot C. On the other hand, 
if the Lya emission from C is purely due to the outflowing 
medium, this doesn't aid in the escape of ionising radiation 
and the LyC photons must be leaking from elsewhere. With- 
out high-resolution LyC imaging, mapping of the Hi appears 
to be the only method by which location of the LyC leakage 
can be addressed. 



5 CONCLUSIONS 

Using HST/ACS we have mapped, calibrated, and anal- 
ysed the Lya emission from nearby luminous blue compact 
galaxy Haro 11. The Lya emission has been compared to: 
Ha; E{B - V) as determined through both Ha/H/3 and UV 
continuum; the kinematic structure of the galaxy from pre- 
vious studies; and the properties of the super star clusters. 
We have used archival X-ray data to map the hot outflow- 
ing gas and to put constraints on the evolutionary state of 
the burst. From SED fltting we have estimated and mapped 
the predicted Lyman-continuum production at 900A, and 
compared this direct measurements. Most notably: 

• Our photometry reproduces spectroscopic fluxes as de- 
termined from the /t/i5 satellite. The total Lya flux is found 
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to be 79.6 x 10 erg s ^ cm ^, corresponding to a lumi- 
nosity of 7.22 X 10*^ erg s"^ 

• The escaping fraction of Lya photons is found to be 
~ 3%. 

• Lyof shows almost no spatial correlation with Ha. The 
Lya morphology shows a central bright emission region 
where Lya photons escape the galaxy directly, surrounded 
by a low surface brightness diffuse halo that results from 
multiple resonant scatterings. 90% of the Lya output is in 
the halo component. Lya/Ha at the values predicted by re- 
combination theory are seen only in the brightest central 
regions. 

• Little correlation is seen between Lya and dust. The 
main central Lya emitting region shows bright emission from 
a region where E{B~V) is significantly greater than regions 
that show only Lya absorption. In the regions of halo emis- 
sion, E{B — V) extends beyond 1 although Lya is resonantly 
scattered and cannot feel the effect of such high extinction. 
This again indicates that dust is not the major regulatory 
factor governing the Lya morphology, which appears to be 
driven more by the Hi distribution and its kinematic struc- 
ture. 

• X-ray observations reveal a diffuse, soft component cen- 
tred on the central star forming knots, and showing hot 
wind-driven gas pushed outside the central starburst region 
as probed by UV bands. The brightest X-ray regions are 
found to be spatially coincident with the regions of highest 
surface brightness in Lya. This indicates that Lya emission 
from central regions may be the result of a perturbed and 
outflowing ISM, accelerated by the release of mechanical en- 
ergy from the SSCs located around knot C. 

• From the super star clusters, we find that peak Lya 
equivalent widths are small at the youngest ages (1 Myr), 
reach a maximum from clusters in the range 2.5-4 Myr, and 
decline again to zero at ages > 8 Myr. This is qualitatively 
consistent with models of Lya emission resulting from out- 
flows in the neutral ISM. 

• Haroll is the only galaxy known to emit Lyman- 
continuum radiation (A < 912A) in the nearby uni- 
verse. From fitting spectral synthesis models we estimate 
the continuum flux at 900A to be /goo.o = 12.3 x 
10~^* erg s"'^ cm~'^ A^^. This corresponds to an escape 
fraction of 9% at 900 A. 
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Table Al. HST/ACS/SBC/FI40LP aperture corrections (in 
magnitudes) derived from TinyTim at radii between 0.10 and 1.0" 
for various power-law continuum slopes. 
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APPENDIX A: SOLAR BLIND CHANNEL 
FI40LP APERTURE CORRECTIONS 

Synthetic PSFs were produced using five different power-law 
continuum slopes {f\ oc A*^) for /3 G [ — 2, —1,0, +1, + 
2 ]. The PHOT task in IRAF is used for aperture photometry. 
Naturally, at large apertures there is little or no variation in 
the correction as a function of (3, although variations exist 
at the 0.014 magnitude level between /3 = —2 and +2. 

This paper has been typeset from a T^JC/ BT^ file prepared 
by the author. 



